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We review recent progress in our theoretical understanding of strongly correlated 
fermion systems in the presence of disorder. Results were obtained by the applica- 
tion of a powerful nonperturbative approach, the Dynamical Mean-Field Theory 
(DMFT), to interacting disordered lattice fermions. In particular, we demon- 
strate that DMFT combined with geometric averaging over disorder can capture 
Anderson localization and Mott insulating phases on the level of one-particle cor- 
relation functions. Results are presented for the ground-state phase diagram of 
the Anderson-Hubbard model at half filling, both in the paramagnetic phase and 
in the presence of antiferromagnetic order. We find a new antiferromagnetic metal 
£^ . which is stabilized by disorder. Possible realizations of these quantum phases with 

• ' ultracold fermions in optical lattices are discussed. 

(N 
O 
O . 

1.1. Introduction 

In non-interacting quantum systems with disorder, e.g., in the presence of randomly 
' distributed impurities, wavefunctions can either be spatially extended or localized. 

j3 ■ Until 1958 is was believed that a localized state corresponds to a bound state of 

an electron at the impurity. By contrast, in his landmark paper of 1958 AndersorP 
predicted that disorder can lead to quite a different type of localized state now re- 
ferred to as "Anderson localized state" . To understand its physical origin it should 
be noted that if a particle is inserted into a disordered system it will start to spread. 
As a consequence the wave is backscattered by the impurities, leading to character- 
istic "weak localization" effectsPSl The multiple scattering of the electronic wave 
can enhance these pcrturbative effects to such a degree that the electrons become 
spatially localized; for reviews seeP^ In this case there is a finite probability for 
an electron to return to the point where it was inserted. If states are extended, this 
probability is zero. So, in contrast to localized states bound at an impurity, Ander- 
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son localized states are confined to a region of space due to coherent backscattering 
from randomly distributed impurities. 

In the thermodynamic limit the excitation spectrum determined from the resol- 
vent of the one-particle system or the one-particle Green functions is very different 
for extended and localized states. The one-particle Green function describing an 
extended state has a branch cut on the real axis, and the spectrum of the Hamilto- 
nian is continuous. By contrast, the Green function for a localized state has discrete 
poles located infinitely close to the real axis, which implies a discrete point spectrum 
of the Hamiltonian. In particular, the point-like spectrum of an Anderson localized 
state is dense. 

In the presence of interactions between the electrons the same classification of 
(approximate) eigenstates may, in principle, be used. Namely, if the one-particle 
Green function of the interacting system has a branch cut at some energies, the 
states at those energies are extended. If the Green function has discrete, separate 
poles the corresponding states are bound states, and if the poles are discrete and lie 
dense the states are Anderson localized. Since one-particle wave functions are not 
defined in a many-body system they cannot be employed to describe the localization 
properties of the system. Instead the reduced one-particle density matrix, or the 
one-particle Green function G(r — r') in position representation, may be employed. 
For localized states these quantities approach zero for |r — r'| — > oo. For extended 
states their amplitude only fluctuates very weakly, i.e., of the order 1/V, where V 
is the volume of the system. 

In the following we are interested in the question how states of many-body 
systems change when the interaction and/or the disorder are varied. In general, the 
very notion of a metal or an insulator is related to the properties of two-particle 
Green functions, e.g., the current- and density-correlation functions. There exist 
different approaches to study the disappearance of a diffusion pole at the metal- 
insulator transition, and correspondingly, the vanishing of the DC conductivity in 
the thermodynamic limit pisEm Q n physical grounds it is very plausible to expect 
that the presence of Anderson localized states with dense, point- like spectrum at the 
Fermi level, discussed above in terms of one-particle Green function, implies zero 
conductivity. Mathematical proofs of this conjecture exist only for specific models 
and in limiting casesP Indeed, it is usually assumed that the presence of Anderson 
localized states at the Fermi level implies the system to be an Anderson insulator, 
at least in the non-interacting case. This is also our line of approach which will be 
reviewed in this article. 

The paper is structured as follows: In Section 2 we review general aspects of the 
interplay between interactions and disorder in lattice fermion systems. In particu- 
lar, we discuss the important question concerning the appropriate average over the 
disorder, and describe the new developments in the field of cold atoms in optical 
lattices which will make it possible in the future to investigate disordered, interact- 
ing lattices fermions with unprecedented control over the parameters. The models 
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of correlated fermions with disorder are introduced in Section 3, followed by an 
introduction into the dynamical mean-field theory (DMFT) (Section 4) and a more 
detailed discussion of arithmetic vs. geometric averaging over the disorder (Sec- 
tion 5). In Section 6 the DMFT self-consistency conditions for disordered systems 
are introduced. After having defined the characteristic quantities which help us to 
identify the different phases of the Anderson- Hubbard Hamiltonian (Section 7) , the 
results for the ground state phase diagram at half filling are reviewed (Section 8). 
In Section 9 the results are summarized. 



1.2. Interplay between interactions and disorder in lattice fermion 
systems 

1.2.1. Interactions vs. disorder 

The properties of solids are strongly influenced by the interaction between the elec- 
trons and the presence of disorder! 4 ^ 8 ! Namely, Coulomb correlations and random- 
ness are both driving forces behind metal-insulator transitions (MITs) which involve 
the localization and derealization of particles. While the electronic repulsion may 
lead to a Mott-Hubbard MITpEGO t ne coherent backscattering of non-interacting 
particles from randomly distributed impurities can cause Anderson localization!^ 
Since electronic interactions and disorder can both (and separately) induce a 
MIT, one might expect their simultaneous presence to be even more effective in 
localizing electrons. However, this is not necessarily so. For example, weak dis- 
order is able to weaken the effect of correlations since it redistributes states into 
the Mott gap and may thus turn an insulator into a (bad) metal. Furthermore, 
short-range interactions lead to a transfer of spectral weight into the Hubbard sub- 
bands whereby the total band-width and thus the critical disorder strength for the 
Anderson MIT increases, implying a reduction of the effective disorder strength. 
Hence the interplay between disorder and interactions leads to subtle many-body 
effects! annmm which pose fundamental challenges for theory and experiment not 
only in condensed matter physics, ^ * 7 * 8 * 18 * 19 ! but most recently also in the field of cold 
atoms in optical lattices pHHUl Indeed, ultracold gases have quickly developed into 
a fascinating new laboratory for quantum many-body physics P * 21 * 28 * ^ A major 
advantage of cold atoms in optical lattices is the high degree of controllability of 
the interaction and the disorder strength, thereby allowing a detailed verification 
of theoretical predictions. The concepts, models, and techniques for their solution 
to be discussed in this paper equally apply to electronic systems and cold fermionic 
atoms in optical lattices. In the following we will therefore refer generally to the 
investigation of "correlated lattice fermion systems" . 
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1.2.1.1. Average over disorder 

In general the theoretical investigation of disordered systems requires the use of 
probability distribution functions (PDFs) for the random quantities of interest. In- 
deed, in physical or statistical problems one is usually interested in "typical" values 
of random quantities which are mathematically given by the most probable value 
of the PDFP However, in many cases the complete PDF is not known, i.e., only 
limited information about the system provided by certain averages (moments or 
cumulants) is available. In this situation it is very important to choose the most 
informative average of a random variable. For example, if the PDF of a random 
variable has a single peak and fast decaying tails the typical value of the random 
quantity is well estimated by its first moment, known as the arithmetic average 
(or arithmetic mean). But there are many examples, e.g., from astronomy, the 
physics of glasses or networks, economy, sociology, biology or geology, where the 
knowledge of the arithmetic average is insufficient since the PDF is so broad that 
its characterization requires infinitely many moments! 35 * 36 ! Such systems are called 
non-self-averaging. One example is Anderson localization: when a disordered sys- 
tem is close to the Anderson MITPmost electronic quantities fluctuate strongly and 
the corresponding PDFs possess long tails which can be described by a log-normal 
distribution!^^ This is well illustrated by the local density of states (LDOS) of 
the disordered system. Most recently it was shown for various lattices in dimen- 
sions d = 2 and 3 that the system-size dependence of the LDOS distribution is an 
unambigous sign of anderson localization, and that the distribution of the LDOS of 
disordered electrons agrees with a log-normal distribution over up to ten orders of 
magnitude.^ Therefore it is not surprising that the arithmetic mean of this ran- 
dom one-particle quantity does not resemble its typical value at all. In particular, 
it is non-critical at the Anderson transitiorPsHUl anc i hence cannot help to detect 
the localization transition. By contrast the geometric mearl 3 - 5 ^ 6 * 46 * 47 * of the LDOS, 
which represents the most probable ("typical") value of a log- normal distribution, 
is the appropriate average in this case. It vanishes at a critical strength of the disor- 
der and hence provides an explicit criterion for Anderson localization in disordered 
systems} 1 * 39 * 42 * even in the presence of interactions! 48 * 49 * 

1.2.1.2. Dynamical mean-field approach to disordered systems 

In general MITs occur at intermediate values of the interaction and/or disorder. 
Theories of MITs driven by interaction and disorder therefore need to be non- 
perturbative. Usually they cannot be solved analytically, and require numerical 
methods or self-consistent approximations. A reliable approximate method for the 
investigation of lattice fermions with a local interaction is provided by dynamical 
mean-field theory (DMFT) JSQHUl where the local single-particle Green function is 
determined self-consistently. If in this approach the effect of local disorder is taken 
into account through the arithmetic mean of the LDOS^l one obtains, in the ab- 
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sence of interactions, the well-known coherent potential approximation (CPA)P^ 
However, the CPA docs not describe the physics of Anderson localization since, as 
discussed above, the arithmetically averaged LDOS is non-critical at the Anderson 
transition.^ To overcome this deficiency Dobrosavljcvic and Kotliai^ formulated a 
variant of the DMFT where the probability distributions (and not only the averages) 
of the local Green functions are determined self-consistcntly ("Statistical DMFT"). 
Employing a slave-boson mean-field theory as impurity solver they investigated the 
disorder-driven MIT for infinitely strong repulsion off half-filling. This statistical 
approach was also employed in other studies of the Hubbard model^ as well as 
in the case of electrons coupled to phonons^ and the Falicov-Kimball modelP^ 
Subsequently, Dobrosavljevic, Pastor, and NikolitpSl incorporated the geometrically 
averaged LDOS into the self-consistency cycle and thereby derived a mean-field 
theory of Anderson localization which reproduces many of the expected features 
of the disorder-driven MIT for non-interacting fermions. This scheme employs 
only one-particle quantities and is therefore easily incorporated into the DMFT for 
disordered electrons in the presence of phononsf^ or Coulomb correlations PSE^D 

1.2.2. Cold atoms in optical lattices: a new realization of disor- 
dered, correlated lattice quantum gases 

During the last few years cold atoms in optical lattices have emerged as a unique 
tool-box for highly controlled investigations of quantum many-body systems. In 
recent years the level of control in applying disordered potentials to ultracold quan- 
tum gases has greatly improved! 23 ^ 24 ! Anderson localization in its pure form has 
been demonstrated by the expansion of weakly interacting Bose-Einstein conden- 
sates in a disordered speckle light field, giving rise to characteristic localized con- 
densate wave functions with exponentially decaying tailsP^SSl The additional in- 
fluence of strong repulsive interactions has been investigated recently in the first 
full experimental realization of the 3d disordered Bose-Hubbard model, by using 
a fine-grained optical speckle field superimposed by an optical latticed In this 
experiment a strong reversible suppression of the condensate fraction due to disor- 
der was observed, indicating the formation of a disorder-induced insulating state. 
Independent experimental evidence was obtained from interacting Rb bosons in 
a quasi-random (bichromatic) optical lattice, where a strong reduction of the Mott 
gap was found and interpreted as possible evidence for a compressible Bose glass 
phased On the theoretical side, low-dimensional quasi-disordered Bose systems 
have been successfully described by DMRG simulations,^ which extended previous 
weak-coupling calculations and found a direct transition from superfluid to Mott 
insulator. Regarding disordered bosons in higher dimensions, the status of theory is 
still more controversial, although significant insight was gained by a new stochastic 
mean-field theory^ which allows for an efficient description of the Bose glass phase 
and has already provided phase diagrams for realistic speckle-type disorder^ such 
as used experimentally.^ Under debate remains the issue of a direct transition 
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between Mott insulator and superfluid, which was claimed to be ruled out in re- 
cent QMC simulations in three spatial dimensions, supported by general heuristic 
arguments.^ Regarding disordered fermions, while no experiments in cold gases 
have been performed yet, theory has significantly advanced in recent years, mostly 
due to progress in the application of DMFT to disordered and inhomogencous sys- 
tems P ^ 9 * 58 * 60 * 66 ! The phase diagram of spin- 1/2 lattice fermions in a random poten- 
tial has now been determined theoretically, both in the paramagnetic phase where 
Mott- and Anderson- insulator compete,^ and in the low-temperature regime where 
antiferromagnetic ordering sets in and a new disorder-induced antiferromagnetic 
metallic phase was foundP In this way, also predictions for single-particle spec- 
tral properties were obtained, which are now becoming accessible experimentally 
via radio frequency spectroscopy measurements of strongly interacting fermionic 
quantum gases in analogy to photoemission spectroscopy of electronic solids. An 
alternative route towards single-particle spectroscopy based on stimulated Raman 
transitions has been discussed theoretically.^ Very recently, also the dynamical 
structure factor of strongly interacting bosons in an optical lattice has been mea- 
sured via two-photon Bragg scattering! 69 * 70 ! These new developments open the door 
towards controlled experimental realization and spectroscopy of strongly interacting 
and disordered fermions in optical lattices. 

1.2.3. Schematic phase diagram 

The Mott-Hubbard MIT is caused by short-range, repulsive interactions in the pure 
system and is characterized by the opening of a gap in the density of states at the 
Fermi level. By contrast, the Anderson MIT is due to the coherent backscattering 
of the quantum particles from randomly distributed impurities in a system with- 
out interactions; at the transition the character of the spectrum at the Fermi level 
changes from a continuous to a dense point spectrum. Already these two limits 
provide great challenges for theoretical investigations. It is an even greater chal- 
lenge to explore the simultaneous presence of interactions and disorder in lattice 
fermions systems. In view of the construction of the dynamical mean-field approach 
employed here, the results which will be presented in the following are expected to 
provide a comprehensive description for systems in spatial dimensions d = 3 and 
larger, i.e., above the limiting dimension d = 2. Two particularly interesting ques- 
tions are whether the metallic phase, which exists at weak enough disorder and/or 
interaction strength, will be reduced or enlarged, and whether the Mott and Ander- 
son insulating phases are separated by a metallic phase. Corresponding schematic 
phase transition lines are shown in Fig. 11.11 It is plausible to assume that both 
MITs can be characterized by a single quantity, namely, the local density of states. 
Although the LDOS is not an order parameter associated with a symmetry breaking 
phase transition, it discriminates between a metal and an insulator which is driven 
by correlations and disorder. 
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Fig. 1.1. Schematic diagram of the possible phases and shapes of phase transition lines in dis- 
ordered, interacting lattice fermion systems. In principle, the interplay between interactions and 
disorder could reduce the metallic regime (dash-dotted line), or enlarge it (full and dashed lines). 
As will be discussed below, investigations within DMFT find that the metallic phase increases 
when interactions and disorder are simultaneously present (full line), and that the two insulating 
phases are connected continuously, i.e., without critical behavior. Insets show the local density of 
states in the absence of disorder or interaction, respectively. 

1.3. Models of correlated fermions with disorder 

Here we study models of correlated fermions on ionic crystals or optical lattices in 
the presence of diagonal (i.e., local) disorder within a tight-binding description. In 
general, these models have the form 

H = *«<r4a c j* + e * n ™ + U Y1 n ^ n ^ ( L1 ) 

ija ia i 

where c\ a and Ci a are the fermionic creation and annihilation operators of the par- 
ticle with spin a = ±1/2 at the lattice site i, n% a = c\ a Ci<, is the particle number 
operator with eigenvalues or 1, and tij a is the probability amplitude for hopping 
between lattice sites i and j. In the Hubbard model Uj^ = iy, i.e. the hopping 
amplitudes are the same for both spin projections. In the Falicov-Kimball model 
tija = tijS a f, i.e. only particles with one spin projection are mobile and the others 
are localized. The second term in (jl.lj) describes the additional external poten- 
tial ej, which breaks the ideal lattice symmetry. For homogeneous systems we set 
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6j = 0, which defines the zero of the energy scale. The third term, a two-body term, 
describes the increase of the energy by U > if two fermions with opposite spins 
occupy the same site. In Eq. only a local part of the Coulomb interaction is 

included and other longer-range terms are neglected for simplicity. Note that this 
approximation is excellent in the case of cold gases in optical lattices, where the 
interaction between neutral atoms is essentially localPSl The disorder affects the 
system through a local term Y) icr ein^, where ti is a random variable drawn from 
a probability distribution function (PDF) V(ei, ...,ejv L ), where Nl is a number of 
lattice sites. Typically we consider uncorrelated, quenched disorder, where 



Each of the P(ei) is the same, normalized PDF for the atomic energies e». The 
quenched disorder means that P(e,) is time independent. In other words, the atomic 
energies are randomly distributed over the lattice and cannot fluctuate in time. 
This type of disorder is different from annealed disorder where the random atomic 
energies have thermal fluctuations. 

In the following we use the continuous box-type PDF 



with Q(x) as the Heaviside step function. The parameter A is therefore a measure 
of the disorder strength. The use of a different continuous, normalized function for 
the PDF would bring about only quantitative changes. 

The Hubbard model and the Falicov-Kimball model defined by (jl.ip are not 
only of interest for solid-state physics, but also in the case of ultracold atoms, where 
specific experimental realizations have been proposedP^l By preparing a mixture of 
bosonic 87 Rb and fermionic 40 K in a 3d optical lattice, Ospelkaus et al. and Giintcr 
et alW^ were able to create - to a first approximation - a version of the Falicov- 
Kimball model where the heavier bosonic species could be slowed down even further 
by using a species-dependent optical lattice and thus become "immobile" while 
the fermionic species remains mobile. Alternatively, if the heavy bosonic species 
could be frozen in a random configuration, this system would allow for a realization 
of the Fermi-Hubbard model with quenched binary onsite disorder. A different 
approach towards quenched randomness in optical lattices was taken by White et 
a?P3 who implemented a fine-grained optical speckle potential superimposed onto 
a 3d optical lattice with interacting bosons and thus realized the bosonic version of 
the Anderson-Hubbard model with continuous disorder. A third alternative 
approach to disordered cold gases is based on bichromatic optical lattices which are 
quasiperiodic, as implemented for the 3d Bose-Hubbard model by Fallani et aZp2l 
who observed a disorder-induced reduction of the Mott excitation gap, similar as 
discussed in the following for the fermionic case. 



V{e 1 ,...,e NL )=\[P{e i ). 



(1.2) 



P(e i ) = - K Q(--\e i \), 



(1.3) 
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The Hamiltonian (jl.ip is not solvable in general. Without disorder, i.e., for 
A = 0, exact solutions on an arbitrary lattice and in arbitrary dimension exist only 
for U = (non-interacting fcrmions), or ij J(T = (fcrmions in the atomic limit). In 
the U — case the solution is obtained via discrete Fourier transform, i.e., 

H = £>- c Lc kCT , (1.4) 

krj 

where = YljCi)^ija- e ~ are f ree fermion dispersion relations. In the 

thermodynamic limit Nl — > oo the spectrum is continuous and eigenstates are 
extended. In the £y ff = limit the lattice sites are uncorrelated and the exact 
partition function has the form Z = Y[ t Zi , where 

Z t = l + 2e^ + er 0u , (1.5) 

where \x denotes the chemical potential within the grand canonical ensemble, and 
(3 = l/ksT is the inverse temperature. In the thermodynamic limit the spectrum 
is point-like and the eigenstates are localized. 

For finite disorder (A ^ 0) an exact solution of the Hamiltonian Ijl.ljl exists 
only for tij a = 0. For a given realization of disorder, i.e., when all values of 
{ei, 62, ejv L } are fixed, the partition function of the model Ijl.ljl is given by 

Z = Uz i = H(l + 2e-^-^ + e-f* u ) . (1.6) 

i i 

As in the atomic limit discussed above {b%j a = 0) the spectrum is point-like in the 
thermodynamic limit and the eigenstates are localized. 

The non-interacting limit (U = 0) of Ijl.ljl with t^ g ^ and disorder A / 
is not exactly solvable. In a seminal paper by Abou-Chacra, Thouless, Anderson^ 
the model Ijl.ljl with U = and tij a = t between nearest neighbor sites was solved 
on the Bethe lattice, which is a tree-like graph without loopsP^EH The solution is 
expressed by the one-particle Green function 

G«( W ) - = 1 — r r t (1.7) 

where the hybridization function 

%M = E ~ — r^r (1-8) 

^ OJ- Ej-Tjjioj) 

describes a resonant coupling of site i with its neighbors. If in the thermodynamic 
limit the imaginary part of rji(z) is finite in some band of energies z, then the states 
with energies z are extended. Otherwise, if the imaginary part of rji(z) is finite at 
discrete energies z such states are localized. For bound states these energies z form 
a point spectrum, and for Anderson localized states the energies z form a dense 
point-like spectrum in the thermodynamic limit. The analysis of the self-consistent 
equations derived for rji(z) by Abou-Chacra, Thouless, Anderson^ showed that, 
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indeed, continuous and dense point spectra are separated by a mobility edge which 
depends on the value of the disorder A. 

In the following we solve the full Hamiltonian (jl.lj) by applying a dynamical 
mean-field approximation to deal with the interaction and then discuss how to cope 
with disorder. 

1.4. Dynamical mean-field theory 

The dynamical mean-field theory (DMFT) started from the following observation!^ 
If the hopping amplitudes are scaled with fractional powers of the space dimension 
d (or the coordination number Z), i.e., t = t* j\[2d = t* j\fZ for nearest neighbour 
hopping on a hypercubic lattice, then in the limit d — > oo (Z — > oo) the self-energy 
£jj(w) in the Dyson equation 

GV^iWn) -1 = G" iCT (iw rl ) _1 - T, ljr7 (iuj n ), (1.9) 

(here in a real-space representation) becomes diagonal 

Y, ija (iuj n ) = Y, icr (iui n ) 5ij, (1-10) 

where lo„ = (2n + l)7r/ /3 are fermionic Matsubara frequencies. In a homogeneous 
system the self-energy is site independent, i.e., T,ij a (iijj n ) = Ti a (iu> n ) 5ij, and is 
only a function of the energy. The DMFT approximation when applied to finite di- 
mensional systems neglects off-diagonal parts of the self-energy. In other words, the 
DMFT takes into account all temporal fluctuations but neglects spatial fluctuations 
between different lattice sites 1S11S21 

Here we apply the DMFT to correlated fermion systems with disorder. Within 
DMFT we map a lattice site onto a single impurity, which is coupled to the dy- 
namical mean-field bath. This coupling is represented by the hybridization function 
r]i a (uj), which is determined self-consistcntly. The mapping is performed for all Nl 
lattice sites. 

The partition function for a particular realization of disorder {ei, €2, £n l } 
is now expressed as a product of the partition functions which arc determined for 
each impurity (representing lattice sites), i.e., 



z =n z '=n cx p(E in [ i 



(i-<H- 7} la (u> n ) - T, la {uj n )] . (1.11) 



The mean- field hybridization function r]i a (u>„) is formally a site- and time- depen- 
dent one-particle potential. In the interaction representation, the unitary time 
evolution due to this potential is described by the local, time-dependent evolution 
operator^EH 

where Ci a (r) evolves according to the atomic part H l ° c of the Hamiltonian (jl.ip in 
imaginary Matsubara time r £ (0, 0), and T T is the time ordering operator. Wc 
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write the partition function (jl.lip as a trace over the operator 

N L 



Z = Z[r ha ] = []Tr [e-^-^^U^]] , (1.13) 



rloc 



where N° c is the local particle number operator. 

Eq. ()1.13|) allows us to determine the local one-particle Green function Gu a {oJn) 
for a given dynamical mean-field 7?io-(w n ). Indeed, the local Green function is ob- 
tained by taking a functional logarithmic derivative of the partition function (|1.13|) 
with respect to r)i a (u) n ), i.e., 

G, CT K) = -^^1. (1.14) 

Then wc find the local Dyson equations 

£ iCT (w„) = iu n + /j,-€i- rj ia (u) n ) - — — - — - , (1-15) 

for each Nl lattice sites. For a single realization of disorder {ei, £2, e« t }, 
Eqs. ([1. 9111.10111. 131 IT. 14) and ll,15| constitute a closed set of equations. A solution 
of this set represents an approximate solution of the Hamiltonian 



1.5. Arithmetic vs. geometric averaging 

A solution of Eqs. (JHH fTTTOl H~T3l H~T4l andHHSJ) is very difficult to obtain in prac- 
tice. For each of the Nl impurities we need to determine the evolution operator 
(|1.12|) exactly Using rigorous methods this can be done only for small Nl- How- 
ever, Eqs. (|1.91 11.101 11131 11.14[ an d I1.15j) should be solved in the thermodynamic 
limit, Nl — > 00. This latter requirement might be overcome by performing finite 
size scaling analysis. But such an analysis requires a large number of lattice sites 
Nl to reliably distinguish Anderson localized states from those belonging to the 
continuum. Here one faces a typical trade-off situation in computational physics. 
The computational problem is greatly reduced when the local interaction in (jl.ljl 
is factorized as in a Hartrec-Fock approximation, whereby genuine correlations are 
eliminated! 15 ! 76 ^ Such approximate treatments can nevertheless provide valuable 
hints about the existence of particular phases. In our invcstigatiorP^HlSl we em- 
ployed the DMFT to include all local correlations as will be discussed in the next 
Section. 

If one could solve the DMFT equations exactly, one would obtain a set of local 
densities of states (LDOS) 

A ia (w) = ——IvaGn a {u>n -> uj + i0+) , (1.16) 

which are random quantities depending on the particular disorder realization 
{ei, £2, £n l }- Usually one needs information about a system that does not 
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depend on a particular disorder realization. Therefore one needs a statistical inter- 
pretation of the solutions of Eqs. (HI IlTOI ITTT51 Ol and EH]). 

When the system is large (cf., Nl — > oo in thermodynamic limit) one usually 
takes the arithmetic average of the LDOS Ai a (u>) over many realizations of the 
disorder, i.e., 

„ n l 

{A ic (u)))= /[[dej P{ei) ^(«;{ eil .., £% }), (1.17) 

J 3=1 

where the dependence on {ei, e2, £n l } is written explicitly. However, such a 
method holds only if the system is self-averaging. This means that samplc-to-samplc 
fluctuations 

D NL {A ia {u)) = (1.18) 

vanish for Nl — > oo, which is equivalent to the Central Limit Theorem for in- 
dependent random variables Ai a (u>). By performing the arithmetic average one 
restores the translational invariance in the description of the disordered system, 
i.e., A CT (a;) ar ith = (Ai a (ui)) is the same for all lattice sites. 

An example of a non-self-averaging system is a disordered system at the Ander- 
son localization transition, or a system whose localization length is smaller than the 
diameter of the sample.^ It implies that during the time evolution, a particle cannot 
explore the full phase space, i.e., cannot probe all possible random distributions. In 
such a case the arithmetic average (| 1 . 1 7[) is inadequate. Here one is faced with the 
question concerning the proper statistical description of such a system. 

The answer was given by Anderson^ One should investigate the full PDF for 
a given physical observable P[Ai a {uj)\ and find its most probable value, the "typ- 
ical" value A CT (a;)typ, for which the PDF P[Ai a (uj)] has a global maximum. The 
typical value of the LDOS, A a (uj) typ , is the same for all lattice sites. By employ- 
ing A a (ui)t yv one restores translational invariance in the description of a disordered 
system. This value will represent typical properties of the system. Using photoe- 
mission spectroscopy one could, in principle, probe the LDOS at a particular lattice 
site and measure its most probable value. We note that if sample-to-sample fluc- 
tuations are small, the typical value A a (uj)typ would coincide with the arithmetic 
average A (T (w) ar ith- On the other hand, in a non-self-averaging system the PDF can 
be strongly asymmetric, with a long tail, in which case the typical value A a (ijj)typ 
would be very different from j4 CT (w) a rith- The arithmetic mean is strongly biased 
by rare fluctuations and hence does not represent the typical property of such a 
system. 

Statistical approaches based on the computation of the probability distribution 
functions would require the inclusion of very many (perhaps infinitely many) im- 
purity sites. This is very hard to achieve in practice, in particular, in correlated 
electron systems discussed here, although there have been recent successful attempts 
in this direction.^ Therefore, one should look for a generalized average which yields 
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the best approximation to the typical value. Among different means the geomet- 
ric mean turns out to be very convenient to describe Anderson localization. The 
geometric mean is defined by 

At M gcom = cxp [(In A ia (w))] , (1.19) 

where (F{ei)) = J Y[ i dei'P(e i )F(ei) is the arithmetic mean of the function F(ei). 
The geometric mean is an approximation to the most probable, typical value of the 
LDOS 

ArMtyp « AxMgcom- (1-20) 

It is easy to see that if P[Ai a (Lu)} is given by a log-normal PDF then A„(uj)t TP = 
A CT (w) g eom holds exactly. It was shown that in the non- interacting case A cr (w) geom 
vanishes at a critical strength of the disorder, hence providing an explicit criterion 
for Anderson localization! 1 ! 39 ! 4 ^ !^ 9 -! We also note that by using the geometrically 
averaged LDOS we restore the translational invariance in our description of a dis- 
ordered system. In addition, as we shall see in the next Section, the restoration 
of translational invariance by averaging allows us to solve the DMFT equations in 
the thermodynamic limit. The problem of finite-size effects is then automatically 
absent. 



1.6. DMFT self-consistency conditions for disordered systems 

According to the spectral theorem the geometrically averaged local Green function 
is given by 

G.K) gcom = / cfc / jHsc ° m . (1.21) 



The DMFT self-consistency condition (|1.14|) is modified now to a translationally 
invariant form 

£ CT (w„) = ico n + 1-1- i] a (uj n ) - 1 . (1.22) 

^tr \yin J geom 

Here we assumed that (e.;) = 0, which holds in particular for the box-shape PDF. 
We also used the translationally invariant hybridization function r] a (uj n ). We can 
now perform a Fourier transform of the lattice Dyson equation (jT79j) and obtain 



G ff K) gcom = / dz- N ° iz) —. , , (1.23) 

J lOJ n - Z + (J, - h a (UJ n ) 

where N (z) is the density of states for a non- interacting and non-disordered lattice 
system. 

Altogether the solution of the DMFT equations for interacting fermions with 
disorder requires the following steps: 

(1) Select (i) Nl values of q from a given PDF P(ei), (ii) an initial hybridization 
function 770- (<*;„), and (iii) an initial self-energy S CT (a; n ); 
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(2) for each solve the impurity problem defined by Eqs. (|1.12l ITTT51 11.141) ; 

(3) determine the LDOS Ai a (u) from the imaginary part of Gu a (uj), and A (7 (i-u) geoTa 
from Eq. ([115]) ; 

(4) employ (|1.21|) to find G <T (w„)g C o m ; 

(5) from Eqs. 11.22"! and ll.23l find a new i] a (Lo n ) and E cr (aj rl ), then go to step (2) until 
convergence is reached. 

It is clear that due to the averaging procedure we restore both translational invari- 
ance and the thermodynamic limit although Nl is finite. Therefore the method is 
superior to other stochastic methods which arc affected by finite size effects. 

In the presence of antiferromagnctic long-range order the self-consistency condi- 
tions are modified. In this case we introduce two sublattices s =A or B, and calculate 
two local Green functions Gu as (uj n ). From this quantity we obtain the geometri- 
cally averaged LDOS A as (u;) seom = exp [(In Ai as (ui))], where Ai as (uj) is given as 
shown in Eq. (|1.16[i . The local Green function is then obtained from the Hilbcrt 
transform (|1.21|) . The local self-energy E ffs (w) is determined from Eq. (|1.22[) . The 
self-consistent DMFT equations are closed by the Hilbert transform of the Green 
function on a bipartite lattice: 



Here s denotes the sublattice opposite to s! 51 | S3 1 

We note that if the geometric mean were replaced by the arithmetic mean one 
would obtain a theory where disorder effects are described only on the level of 
the CPA, which cannot detect Anderson localization. It should also be pointed 



normalized to unity. This means that A a (uj) geom only describes the extended states 
of the continuum part of the spectrum. Localized states, which have a dense point 
spectrum, are not included in the DMFT with geometric average. Therefore, this 
approach cannot describe the properties of the Anderson-insulator phase. 

The accuracy of the DMFT approach with geometric average over disorder was 
checked against numerically exact results obtained for non-interacting fermions on 
a cubic lattice.SSEU The critical disorder strengths at which Anderson localization 
occurs were found to agree within a factor of two^or better PS However, there exists 
a discrepancy regarding the shape of the mobility edge, which shows a pronounced 
reentrant behavior for non-interacting particles with box-type PDF of the disorder. 
This feature is not reproduced by our approach.^ On the other hand, the reentrant 
behavior is a non-universal feature. Namely, it is much less pronounced in the case 
of a Gaussian PDF for disorder, and does not occur at all for a Lorentzian PDFpSl 

It should be pointed out that the DMFT-based self-consistent approach to in- 
teracting lattice fermions with disorder discussed here, is not related to the self- 
consistent theory of Anderson localization by Vollhardt and Wolfle^ and its gen- 
eralizations! 80 * 81 ! Namely, the latter theory determines the frequency dependent 




(1.24) 



out that in the presence of disorder the LDOS represented by A a {uS) 



gcom 



is not 
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diffusion coefficient D(u>) from arithmeticaily averaged two-particle correlation func- 
tions by considering diffuson and cooperon diagrams. The approach reviewed here 
does not make use of these coherent back-scattering contributions, but computes a 
one-particle correlation function, the LDOS, and thereby extracts information on 
Anderson localization. The fact that the DMFT is based on a local approximation 
through the limit of large spatial dimensions does not necessarily imply that back- 
scattering contributions arc entirely absent in this approach. Indeed, contributions 
due to back-scattering are implicitly contained in the hybridization function, which 
describes the diffusion of one-particle excitations away from and back to a given 
lattice siteP^H Quite generally the relation between theoretical approaches based 
on one-particle and two-particle correlation functions, respectively, and their results 
for the critical disorder strength for Anderson localization, is still not sufficiently 
understood and will continue to be an important topic for future research. Perhaps 
the limit of high lattice dimensions will serve as a useful starting point .HH 82 I 83 I 

1.7. Identification of different phases 

To characterize the ground state of the Hamiltonian the following quantities 
are computed: 

(1) the LDOS A as (u}) gcora for a given sublattice s and spin direction er; 

(2) the total DOS for a given sublattice s at the Fermi level (uj = 0) with 

A" S (0)goom = J2a Ars(w = 0) gcO m; 

(3) the staggered magnetization "i^ ™ 1 = |«|™ m — n|^ m |, where nf° m = 
j ^ dujA as ((jj) gcora is the local particle density on sublattice 

For comparison we determine these quantities also with the arithmetic average. 

The possible phases of the Anderson-Hubbard model can then be classified as 
follows: The systems is a 

• paramagnetic metal if A r | eom (0) ^ and m^" 1 = 0; 

• AF metal if Af oom (0) + and mf£ m ^ 0; 

• AF insulator if Af oora (0) = and mf° m ^ but Af cora (w) ^ for some u ^ 
(in fact, the last condition is already implied by rn s £y la 7^ 0); 

• paramagnetic Anderson-Mott insulator if A| oom (w) = for all uj. 

Note, that we use the term "metal" also for neutral fermionic atoms if they fulfil 
the above conditions. 

1.8. Ground state phase diagram of interacting, disordered lattice 
fermion systems at half-filling 

We now apply the formalism discussed above to the Anderson-Hubbard model at 
half-filling and compare the ground state properties in the paramagnetic and mag- 
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Fig. 1.2. Non-magnetic ground state phase diagram of the Anderson-Hubbard model at half- 
filling as calculated by DMFT with the typical local density of states; after Ref. 58. 
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In the following we choose a model DOS, No(e) = 2v 'D 2 — e 2 /ttD 2 , with band- 
width W = 2D, and set W = 1. For this DOS and for a bipartite lattice the local 
Green function and the hybridization function are connected by the simple algebraic 
relation 77 ffs (w) gcom = D 2 G txS (ui) geom / '4.ED 

The DMFT equations are solved at zero temperature by the numerical renormal- 
ization group technique,^ which allows us to calculate the geometric or arithmetic 
average of the local DOS in each iteration loop. 



1.8.1. Paramagnetic phase diagram 

The ground state phase diagram of the Anderson-Hubbard model at half-filling ob- 
tained within the DMFT approach discussed above is shown in Fig. II. Two 
different phase transitions are found to take place: a Mott-Hubbard MIT for weak 
disorder A, and an Anderson MIT for weak interaction U . The correlated, disor- 
dered metal is surrounded by two different insulating phases whose properties, as 
well as the transitions between them, will now be discussed. In this section the spin 
index a is omitted since all quantities are spin independent. 

(i) Disordered, metallic phase:The correlated, disordered metal is characterized 
by a non-zero value of the spectral density at the Fermi level, A(u> = 0) geO m 7^ 0. 
In the absence of disorder DMFT predicts this quantity to be given by the bare 
DOS Nq({)), which is a consequence of the Luttinger theorem. This means that 
Landau quasiparticles are well-defined at the Fermi level. The situation changes 
completely when disorder is introduced since a subtle competition between disorder 
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A U 



Fig. 1.3. Left panel: Local density of states (LDOS) as a function of disorder A for various 
values of the interaction U . Solid (dashed) curves correspond to the geometrically (arithmetically) 
averaged LDOS. Right panel: Geometrically averaged LDOS as a function of interaction U for 
different disorder strengths A. Solid (dashed) curves with closed (open) symbols are obtained 
with an initial metallic (insulating) hybridization function. Triangles: A = 1; dots: A = 2.5. Left 
inset: LDOS with Mott gap at U = 3 for different disorder strengths A. Right inset: Integrated 
LDOS N gcoIn as a function of A at U = 3; after Ref. 58. 

and electron interaction arises. 

Increasing the disorder strength at fixed U reduces A(0) gcom and thereby de- 
creases the metallicity as shown in the upper panel left of Fig. 11.31 The opposite 
behavior is found when the interaction is increased at fixed A (see right panel of 
Fig. II .31 for A = 1), i.e., in this case the metallicity improves. In the strongly inter- 
acting metallic regime the value of A(0) goom is restored, reaching again its maximal 
value iVo(0). This implies that in the metallic phase sufficiently strong interactions 
protect the quasiparticles from decaying by impurity scattering. For weak disorder 
this interaction effect is almost independent of how the LDOS is averaged. 

(ii) Mott- Hubbard MIT: For weak to intermediate disorder strength there is a 
sharp transition at a critical value of U between a correlated metal and a gapped 
Mott insulator. Two transition lines are found depending on whether the MIT is 
approached from the metallic side [A*| ff ([7), full dots in Fig. 11.2] or from the 
insulating side [A^[ H (U), open dots in Fig. 11.2] , The hysteresis is clearly seen in 
right panel of Fig. Ofor A = 1. The curves A*( H (U) and A^ H (U) in Fig. Oare 
seen to have positive slope. This is a consequence of the disorder-induced increase 
of spectral weight at the Fermi level which in turn requires a stronger interaction 
to open the correlation gap. In the Mott insulating phase close to the hysteretic 
region an increase of disorder will therefore drive the system back into the metallic 
phase. The corresponding abrupt rise of A(0) geonl is clearly seen in the left lower 
panel of Fig. 11.31 In this case the disorder protects the metal from becoming a 
Mott insulator. 

Around A w 1.8 the curves A^( H (U) and A^ H (U) terminate at a single critical 
point, cf. Fig. 11.21 For stronger disorder (A > 1.8) there appears to be a smooth 
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crossover rather than a sharp transition from the metal to the insulator. This is 
illustrated by the U dependence of A(0) gC om shown in right panel of Fig. 11.31 for 
A = 2.5. In this parameter regime the Luttinger theorem is not obeyed for any U. 
In the crossover regime, marked by the hatched area in Fig. I1.2[ A(0) goom vanishes 
gradually, so that the metallic and insulating phases can no longer be distinguished 
rigorously 

Qualitatively, we find that the Mott-Hubbard MIT and the crossover region 
do not depend much on the choice of the average of the LDOS.^ We also note 
the similarity between the Mott-Hubbard MIT scenario discussed here for disor- 
dered systems and that for a system without disorder at finite temperatures ISUSU 
especially the presence of a coexistence region with hysteresis. However, while in 
the non-disordered case the interaction needed to trigger the Mott-Hubbard MIT 
decreases with increasing temperature, the opposite holds in the disordered case. 

(Hi) Anderson MIT: The metallic phase and the crossover regime are found to lie 
next to an Anderson insulator phase where the LDOS of the extended states vanishes 
completely (see Fig. II .2|) . The critical disorder strength Af(U) corresponding to 
the Anderson MIT is a non-monotonous function of the interaction: it increases 
in the metallic regime and decreases in the crossover regime. Where A*(U) has 
a positive slope an increase of the interaction turns the Anderson insulator into a 
correlated metal. This is illustrated in Fig. 11.31 for A = 2.5: at U/W ps 0.7 a 
transition from a localized to a metallic phase occurs, i.e., the spectral weight at 
the Fermi level becomes finite. In this case the electronic correlations inhibit the 
localization of quasiparticles by scattering at the impurities. 

Fig. II .31 shows that the Anderson MIT is a continuous transition. In the critical 
regime ^(0)gcom ~ [A^([/) - Af for U = const. In the crossover regime a critical 
exponent /3 = 1 is found (see the case U = 1.75 in lower panel of Fig. II. 3[) : elsewhere 
(5 7^ 1. However, we cannot rule out a very narrow critical regime with f3 = 1 since 
it is difficult to determine /? with high accuracy. It should be stressed that an 
Anderson transition with vanishing A(0) gcom at finite A = A^{U) can only be 
detected within DMFT when the geometrically averaged LDOS is used (solid lines 
in Fig. II. 3[) . Indeed, using arithmetic averaging one finds a nonvanishing LDOS at 
any finite A (dashed lines in Fig. II. 3[) . 

(iv) Mott and Anderson insulators: The Mott insulator (with a correlation gap) 
is rigorously defined only in the absence of disorder (A = 0), and the gapless 
Anderson insulator only for non-interacting systems (U = 0) and A > A^(0). For 
finite interactions and disorder this distinction can no longer be made. On the other 
hand, as long as the LDOS shows the characteristic Hubbard subbands (left inset in 
Fig. II. 3|) one may refer to a disordered Mott insulator. With increasing disorder A 
the spectral weight of the Hubbard subbands vanishes (right inset in Fig. I1.3[) and 
the system becomes a correlated Anderson insulator. The boundary between these 
two types of insulators is marked by a dashed line in Fig. 11.21 The results obtained 
here within DMFT show that the paramagnetic Mott and Anderson insulators are 
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continuously connected. Hence, by changing U and A it is possible to move from 
one insulating state to another one without crossing a metallic phase. 

1.8.2. Magnetic phase diagram 

At half filling and in the absence of frustration effects interacting fermions order an- 
tifcrromagnctically. This raises several basic questions: (i) How is a non-interacting, 
Anderson localized system at half filling influenced by a local interaction between 
the particles? (ii) How does an antiferromagnetic (AF) insulator at half filling re- 
spond to disorder which in the absence of interactions would lead to an Anderson 
localized state? (iii) Do Slater and Heisenberg antiferromagnets behave differently 
in the presence of disorder? Here we provide answers to these questions by calculat- 
ing the zero temperature, magnetic phase diagram of the disordered Hubbard model 
at half filling using DMFT together with a geometric average over the disorder and 
allowing for a spin-dependence of the 

The ground state phase diagram of the Anderson-Hubbard model (jTTTJ) obtained 
by the above classification is shown in Fig. 11.41 Depending on whether the inter- 
action U is weak or strong the response of the system to disorder is found to be 
very different. In particular, at strong interactions, U/W > 1, there exist only two 
phases, an AF insulating phase at weak disorder, A/W < 2.5, and a paramagnetic 
Andcrson-Mott insulator at strong disorder, A/W > 2.5. The transition between 
these two phases is continuous. Namely, the local DOS and the staggered mag- 
netization both decrease gradually as the disorder A increases and vanish at their 
mutual boundary (lower panel of Fig. I1.5p . By contrast, the phase diagram for 
weak interactions, U/W < 1, has a much richer structure (Fig. II. 4| ). In particular, 
for weak disorder a paramagnetic metallic phase is stable. It is separated from the 
AF insulating phase at large U by a narrow region of AF metallic phase. The AF 
metallic phase is long-range ordered, but there is no gap since the disorder leads to 
a redistribution of spectral weight 

To better understand the nature of the AF phases in the phase diagram we take a 
look at the staggered magnetization rn% F - The dependence of m^ ™ 1 on U is shown 
in the upper panel of Fig. ll.5l for several values of the disorder A. In contrast to the 
non-disordered case a finite interaction strength U > U C (A) is needed to stabilize 
the AF long-range order when disorder is present. The staggered magnetization 
saturates at large U for both averages; the maximal values depend on the disorder 
strength. In the lower panel of Fig. 11.51 the dependence of to^f on the disorder A 
is shown for different interactions U. Only for small U do the two averages yield 
approximately the same results. 

Another useful quantity is the polarization P£ F = m\ F /I a : where I a = 
f-x Scts / 9 "s( tlJ )^ w /2 is the total spectral weight of p" s (uj). It allows one to inves- 
tigate the contribution of the point-like spectrum of the Anderson localized states 
to the magnetization. This provides important information about the spectrum 
since with increasing disorder more and more one-particle states of the many-body 
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Fig. 1.4. Magnetic ground state phase diagram of the Anderson-Hubbard model at half-filling as 
calculated by DMFT with a spin resolved local DOS (see text); PM: paramagnetic, AF: antifer- 
romagnetic; after Ref. 60. 

system are transferred from the continuous to the point-like spectrum. For weak 
interactions (U = 0.5) the decrease of the polarization with increasing disorder A 
obtained with geometric or arithmetic averaging is the same (see inset in Fig. 1 1 . 5j) . 
Since within arithmetic averaging all states are extended, the decrease of 
(which is also the same for the two averages in the limit of weak interactions, see 
lower panel of Fig. I1.5[) must be attributed to disorder effects involving only the 
continuous spectrum. At larger U the polarization is constant up to the transition 
from the AF insulator to the paramagnetic Anderson-Mott insulator. In the lat- 
ter phase the polarization is undefined, because the continuous spectrum does not 
contribute to J^ c F ° m . 

In the absence of disorder the AF insulating phase has a small ("Slater") gap 
at U/W < 1 and a large ("Heisenberg") gap at U/W > 1. These limits can be 
described by perturbation expansions in U and 1/U around the symmetry broken 
state of the Hubbard and the corresponding Heisenberg model, respectively. In 
agreement with earlier studies^ our results for toaf (upper panel of Fig. [TT5"]) show 
that there is no sharp transition between these limits, even when disorder is present. 
This may be attributed to the fact that both limits are described by the same order 
parameter. However, the phase diagram (Fig. 11.4)1 shows that the two limits can be 
distinguished by their overall response to disorder. Namely, the reentrance of the 
AF metallic phase at A/W > 1 occurs only within the Slater AF insulating phase. 

The magnetic structure of the Anderson-Mott insulator cannot be determined 
by the method used here since it describes only the continuous part of the spectra 
and not the point spectrum. However, only the paramagnetic solution should be 
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Disorder, A 

Fig. 1.5. Upper panel: Staggered magnetization rn^p™ as a function of interaction U. Lower 
panel: "i^P' a = geom/arith, as a function of disorder A (full lines: arithmetic average, dashed 
lines: geometric average). Inset: Polarization P^ F as a function of disorder!^ Reprinted with 
permission from Phys. Rev. Lett. 102, 146403 (2009). © American Physical Society. 

expected to be stable because the kinetic exchange interaction responsible for the 
formation of the AF metal is suppressed by the disorder. This does not exclude the 
possibility of Griffiths phase-like AF domains! 89 * 90 ! 

It is interesting to note that even the DMFT with an arithmetic average finds 
a disordered AF metalP^^ However, the arithmetically averaged local DOS in- 
correctly predicts both the paramagnetic metal and the AF metal to remain stable 
for arbitrarily strong disorder. Only a computational method which is sensitive to 
Anderson localization, such as the DMFT with geometrically averaged local DOS 
employed here, is able to detect the suppression of the metallic phase for A/W > 1.5 
and the appearance of the paramagnetic Andcrson-Mott insulator at large disorder 
A already on the one-particle level. 

1.9. Summary 

In this article we reviewed the properties of low-temperature quantum phases of 
strongly correlated, disordered lattice fermion systems with application to corre- 
lated electronic systems and ultracold fermions in optical lattices. We discussed the 
Anderson-Hubbard model and a comprehensive nonperturbative theoretical method 
for its solution, the dynamical mean-field theory combined with geometrical aver- 
aging over disorder. This approach provides a unified description of Anderson- and 
Mott-localization in terms of one-particle correlation functions. 

We presented low-temperature quantum phase diagrams for the Anderson- 
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Hubbard model at half filling, both in the paramagnetic and the antiferromag- 
netic phase. In the paramagnet we observed re-entrant metal-insulator transitions 
induced by disorder and interaction, where the corresponding Anderson- and Mott- 
insulating phases are continuously connected. In the presence of antiferromagnetic 
order, a new antiferromagnetic metallic phase was found, which is stabilized by the 
interplay between interaction and disorder. 

It is expected that these new quantum states will be observable by using ul- 
tracold fermions in optical lattices where disorder and interactions are easily tun- 
able in a wide range. While current experimental temperatures are still above 
those required for observing quantum antiferromagnetism, the paramagnetic Mott- 
Anderson insulator should be easily accessible within current setups. 

Even after several decades of research into the complex properties of disor- 
dered, interacting quantum many-body systems many fundamental problems are 
still unsolved. Future investigations of the existing open questions, and of the new 
questions which are bound to arise, are therefore expected to provide fascinating 
new insights. 
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